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Abstract 

The proton momentum and density distributions of closed shell nuclei are calcu- 
lated with a model treating short-range correlations up to first order in the cluster 
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expansion. The validity of the model is verified by comparing the results obtained 
using purely scalar correlations with those produced by finite nuclei Fermi Hy- 
pernetted Chain calculations. State dependent correlations are used to calculate 
momentum and density distributions of ^^C, ^^O, ^°Ca, and ^^Ca, and the effects 
of their tensor components are studied. 



PACS number(s): 21.10.Ft, 21.60.-n 
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1 Introduction 

The high precision data produced by modern electron scattering experiments have im- 
posed severe constraints on the vahdity of the models and the theories aiming to describe 
the nuclear properties. Since the beginning of the eighties, elastic scattering experiments 
|l|], measuring charge distributions, and knock-out experiments [0, measuring spectral 
functions, have pointed out the difficulties of the mean-field model in the description of 
the ground state of atomic nuclei. 

The effects not considered in the mean-field model have been generically named cor- 
relations. One usually distinguishes between short- and long-range correlations. The 
former ones are those acting at short inter particle distances modifying the mean field 
single particle wave functions to take care of the hard core part of the nuclear potential. 
The long-range correlations are instead acting on the full system and are produced by 
collective phenomena like sound waves or surface vibrations. 

The Correlated Basis Function (CBF) theory, based on Jastrow's approach is 
particularly suitable for the study of the short-range correlations, since these correlations 
are explicitly included in the definition of the many-body wave function. Unfortunately 
the technical complexity of this theory has limited its application to few-body systems 
and to infinite nuclear matter |^. At present there are only few example of CBF 
calculations in nuclei heavier than "^He. The CBF theory has been used in Variational 
Monte Carlo calculations of the ^^O nucleus [0. Only recently the Fermi Hypernetted 
Chain Theory (FHNC) has been extended to describe the ground state properties of 
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medium- heavy doubly-closed shell nuclei ^ |T0|- Unfortunately these calculations are 
still limited to the use of central interactions and scalar correlations. 

Because of these difficulties, the study of short range correlations in medium-heavy 
nuclei has been done using simplified models. Nuclear models considering short-range 
correlations have been used to analyse elastic electron scattering data already at the end 



of the sixties [|TT1. These models are based on cluster expansion like the CBF theory, but 



they retain only those terms containing a single correlation line [T^ . The truncation of the 
expansion simplifies the calculations and it is done in a way to conserve the normalization 
of the density distribution. More recently, these simplified models have been used to 



investigate also the nucleon momentum distributions [0, 0, |T5 . 

The aim of the present article is twofold: first we would like to discuss the validity of 
these nuclear models, and secondly we would like to investigate the effects of the state 
dependent part of the correlation on density and momentum distributions. 

For these purposes we extend the model developed in Ref. |TB[ to calculate both 



density and momentum distributions for doubly closed shell nuclei and we compare our 



results with those of finite nuclei FHNC calculations [T^. After this test, using a nuclear 



matter correlation function [|17|, we study the effects produced by the state dependent 
terms of the correlation on the proton density and the momentum distributions of the 
12Q 160 4°Ca and ^^Ca nuclei. 
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2 The model 

The basic quantity of interest for our calculations is the One-Body Density Matrix 
(OBDM) defined as: 

A r 

p(ri, r'l) = j^J d?r2 d^n ... rfV^^'^(ri, rs, rg, r^)^(r;, ra, n, ta), (1) 

where is the nuclear ground state wave function and A/" =< >. In Eq. (|l|) a sum 
on the spin and isospin components of all the particles, particle 1 included, is understood. 
In our model the protons and the neutrons are separately treated, therefore: 

p(ri,r;)=//(ri,r;)+p"(ri,r;) (2) 

where the protons and neutrons OBDMs are obtained inserting in Eq. (|l]) the operators 
selecting the protons, 

Q(l) = ^(1 + ^3(1)), (3) 

and the neutrons, 1 — Q{1). 

The density distribution is the diagonal part of the OBDM (ri = r[) while the mo- 
mentum distribution is its Fourier transform: 

n(k) = jj d?n J d\[p{ri, r[)e'''<'''-<\ (4) 

The OBDM is normalized as: 

J dVi J rfV;p(ri, r[)6{r^ - r[) = A (5) 

and therefore: 

'd^kn{k) = {2'Kf. (6) 



Our model is based upon the Jastrow ansatz on the nuclear wave function 0]: 

^{1,...,A) = F{1,...,A)^1,...,A) (7) 

where F is an A-body correlation operator, and $ is a Slater determinant built on a basis 
of single particle (s.p.) wave functions generated by a one-body hamiltonian. 

In the present work we shall use a spherical mean-field potential with a spin-orbit 
interaction, therefore the quantum numbers characterizing the s.p. wave functions are 
the principal quantum number n, the orbital angular momentum I, the total angular 
momentum j, its third component m and the 3rd isospin component t. 

A usual choice of the expression of the correlation operator is based upon an extension 
of the original Jastrow's ansatz to a state dependent form: 

A M 

Fil,...,A) = S n [Y.fnin,)0-^, (8) 

j>i=l n=l 

where Vij = |rj — rj|, and we have indicated with S a symmetrizer operator, with M the 
maximum number of the correlation channels and with O^^j the operators characterizing 
the various channels. In the present work we use correlations with 6 active channels 
defined as: Oj^ = 1, Of. = a{t) ■ ^t(j), = t(^) ■ t(j) Of^ = ■ aij^i) ■ t(j), 
Ofj = Sij, Ofj = SijT{i) ■ T{j), where Sij is the tensor operator: 

S,=3Wi^llW^ (9) 

The OBDM is calculated with cluster expansion techniques applied to finite systems 
p!8| . We express the correlation function as: 

A M 

F(1,...,A) = 1 + S n [Y.hn{r.,)0?j\, (10) 

j>i=l n=l 
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and perform the expansion in terms of hn- 

The denominator in Eq. cancels all the unlinked diagrams of the numerator, i.e. 
all the diagrams which have not a direct link, either statistical or dynamical one, with 
the points r and r'. 

The basic hypothesis of the model enters now, because of all the infinite set of linked 
diagrams obtained within the cluster expansion, we retain only those containing corre- 
lations lines up to the second order in the correlation function hn. The set of diagrams 
considered is shown in Fig.l. This is the lowest order set of correlated diagrams conserving 
the normalization of the OBDM, Eq. (|^). 

It is possible to describe this limited set of diagrams in terms of the correlation function 
/„. In this new representation the set of diagrams to be calculated is shown in Fig.2. 

The OBDM calculated with the diagrams of Fig.2 can be expressed as: 

p?(ri,r;) = pl{r^,v[) + A{vur[)-Bir^,r[)-Civur[) + D{v^,r[) 
= Po(ri> r'J + po(ri, J d^r2 H{ri, r[, r2)po(r2, 
-/*V,p.(r,..,) HMM,y,) 
-/*V./d'r,p„(r,.r.)p„(r..r;)p„(r3.r3) i/(r..r..r3) 
+ /<;S/<i'r3P„(r„r.)po(r.,r>„(r3,r;) //(r.r.r,). (11) 

where, like in Eq. (|^), the sum on spin and isospin components is understood. 

In the above expression we have used the uncorrelated OBDM defined in terms of the 
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s.p. wave functions as: 

Po(ri,r2)= E CVm(ri)0n«,Jr2). (12) 

nljmt 

In the diagrams of Fig.2 the po(ri, is represented by an oriented hne. The dashed hnes 
represent the dynamical correlations expressed in Eq. ( |Tl| ) by the coefficients H{ri, r2, r^) 
defined as: 

Hir,,V2,r,) = (E/p(n3)0?3)Q(l)(E/,(r23)0?3), (13) 

p=l q=l 

with g(l) defined by Eq. (|). 

It is easy to see that the set of considered diagrams conserves the density normalization. 
Because of the property of the uncorrelated density 

J (frj po(ri, rj)po(ri, ^k) = Po(ri, rfe), (14) 

where we have also summed on the spin and isospin coordinates of the particle j, we 
obtain the following relations among the diagrams of Fig.2: 

J d^n J d^r[ A(ri, r[) 5(ri - r'^) = J d^n J d^r[ C(ri, r[) 5(ri - r[) (15) 

J d^n J d^r[ Bin, r'l) 6{ri - r[) = J d\i J d\[ D{ri, r[) 6{ri - r[) (16) 
and therefore: 

I d'n I d'r[ p?(ri, v[) 5(ri - r'^) = J dh, j dh[ pg(ri, r',) 5{v, - v[) = Z (17) 

In the evaluation of the spin and isospin traces of the correlation kernels of Eq. (|T3p, 
we found convenient to separate the various terms with respect to their isospin dependence 
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and we obtain four terms: 



if(ri,r;,r2) = 

figi + 3/2^/2 + 6/5(/5 + {fi92 + f29i - 2/292 + 2/5^5) o-(l) ■ a{2) + 
{fi95 + h9i + f29b + 1^92 - 2/55(5) S'12] (5(1) + 

fm + 3/4^?2 + 6/6(?5 + {h92 + h9i - 2U92 + 2/6^75)^(1) ■ t(2) + 
(/3^75 + f69i + h95 + h92 - 2hg,)S,2]T{l) ■ t{2)Q{1) + 

'h9z + 3/2(74 + 6/5^76 + [h9A + /2^73 - 2/2(74 + 2/556)«t(1) ■ ct{2) + 
(/i^76 + /5^73 + f29& + /5^74 - 2/5(76) ^12] Q(1)t(1) ■ r(2) + 

/sS-s + 3/45-4 + 6/656 + (/3fl'4 + /4fl'3 " 2/454 + 2^56) cr(l) ■ cr(2) + 
(/3^76 + /6^73 + /4(76 + h9A ' 2/656) ^12] t(1) ■ t(2)Q(1)t(1) ■ t(2) 



1^ 



In the above expression to simplify, the writing, we used the symbol /„ for the correlation 
functions depending from ri and r2, and 5„ for those depending from r'^^ and r2. Clearly, 
in the case ri = r'^ we have /„ = 9n- 



The calculation continues inserting the above expression in Eq. ([lID and evaluating 



the spin and isospin traces. The traces of the isospin dependent part of each term are 
given in Tab.l. We have explicitly verified that the spin traces of all the terms depending 
from the tensor operator are zero. This is due to the definition of the tensor operator, 
Eq. (^, to the spherical symmetry of the problem and to the saturation of the spin of all 
the single particle wave functions. 

In order to express the final result we used the quantity po^'''^'*(ri, r2) defined by the 
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relation: 

Po(ri, r^) = E E Po''"'\^u r2)xli(l)x.2(2)xl(l)x*(2), (19) 

t=p,n sls2 

where we have indicated with x the spin and isospin wave functions. Together with the 
above definition we define as weU: 

Po(ri,r2) = EPr*(ri,r2) (20) 

s 

p*(rO = p*(ri,ri). (21) 
Using the quantities defined above we express the four terms composing p^, Eq. (|TTp, 

as: 

A{r,, r[) = pg(ri, r[) J rfV2[pg(r2)Gl(ri, r',, r^) + p^oi^2)G2{r„ r[, r,)] 

+ 4po"(ri,r;)|dV2pg(r2)G3(ri,r;,r2) (22) 

i?(ri, r[) = f £r,{ J] p^^'^^'^(ri, r,) p^'"^''^ (r,, v[) G4(ri, r',, r^) 

sls2 

+ 2p[;(ri, r2)pg(r2, r'^) G5(ri, r'^, rs) 

+ E 2po^''^'^(ri, r2)po^'"^'"(r2, r[) G6{ri, r[, rs) 

sls2 

+ 4pg(ri,r2)p^(r2,r;) G7{r,,r[,r2) 

+ Y: 2po''^'''^(ri, rM''''''ir2, r[) G8{t„ r[, r,) 

sls2 

+ 4po"(ri,r2)pg(r2,r'i) G'9(ri,r;,r2)} (23) 
C(ri, v[) = f dh, f d'rs E Po''^''"(ri, r2)pf '^^'^(r2, r',) 

sls2 
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pg(r3)Gl(r2, r2, r,) + p^(r3) (G2(r2, r^, r,) + 4G'3(r2, r^, r,))] (24) 



si s2 s3 

[po^''^'^(r2, r3)G'4(r2, r2, rs) + 2po^'''^'"(r2, rg) (G6(r2, r2, rg) + G'8(r2, r2, rg)) 

I sl,s2,p/ \ s2.sl.p/ I \ 

sls2 

2 [pg(r2, r3)G'5(r2, r2, rg) + 2p^(r2, rg) (G'7(r2, v^, rg) + G9(r2, r2, rg))] } (25) 



The explicit expression of the terms G containing the correlation functions are given in 
the Appendix. 

The calculation of the momentum and density distributions of doubly-magic nuclei has 
been done by exploiting the spherical symmetry of the problem to perform an analytical 
integration on the angular variables. For this purpose we have expressed the s.p. wave 
functions in spherical coordinates and we have performed a multipole expansion of the 
functions G containing the correlations. In the appendix we present the basic points of 
the calculation and we give the final expressions used to evaluate the diagrams A, B, G 



and D. The details of the calculation are presented in Ref. [19 



3 Results 

In our calculations the s.p. wave functions have been generated by a spherical Woods- 
Saxon well of the form: 
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where is the pion mass. 

For each nucleus considered, we have taken from the hterature the parameters of 
the potential (see Tab. 2). These parameters have been fixed to reproduce the s.p. energies 
around the Fermi surface and the root mean squared charge radii. 

We tested the validity of our model by comparing our results with those of FHNC 
calculations performed with the same inputs. We used two kinds of correlation functions. 
A first one has a gaussian functional dependence from the inter particle distance. The 
second one has been obtained in Ref. by solving the variational equations for the 
mean value of the hamiltonian up to the second order (Euler correlations). 

The correlations used for the ^^C and '^'^Ca nuclei are shown in Fig. 3. The dashed lines 
represent the gaussian correlations, while the full lines show the Euler correlations. These 



correlations have been fixed in Ref. |10| by minimizing the nuclear binding energies for 



the Afnan and Tang S3 interaction [|2T|. The Euler correlations show an overshooting of 
the asymptotic value in the region between 1 and 2 fm. 

In Fig. 4 we compare the momentum and density distributions calculated within the 
FHNC theory [jlOl (dashed hues) and those obtained with the present model (full lines). 
The agreement between the two calculations is quite good, independently from the cor- 
relation function used. 

Confident of the validity of our model we use the nuclear matter state dependent 



correlation function of Ref. to calculate momentum and density distributions in C, 



^^O, ^oCa and ^^Ca. 



We show in Fig. 5 the various terms of the correlation as a function of the inter particle 
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distance. The dominant term of the correlation is the scalar channel. The tensor term 
(/s) is extremely small and the tensor-isospin term (/g) is peaked at an inter particle 
distance of about 1 fm. 

In Figs. 6 and 7 we show the proton density and momentum distributions for the 
four nuclei considered calculated by switching on and off the various channels of the 
correlation function. The dotted lines represent the uncorrelated results. The results of 
the calculations performed with the correlation containing all the six channels are shown 
by the full lines. The dashed lines have been obtained using correlations with only the 
scalar channel active, while the dashed-dotted lines with correlations without the two 
tensor channels. 

The results of Fig. 6 show the same behaviour for each nucleus considered. In the 
nuclear interior, the correlated density distributions are smaller than the uncorrelated 
ones. We notice that the curves obtained with purely scalar correlations are the more 
distant ones from the uncorrelated results. The inclusion of the other correlation channels 
reduces these differences. 

The momentum distributions (Fig. 7) show a well known behaviour. Correlated and 
uncorrelated results are practically the same up to momentum values of about 1.8-2.0 
fm~^, but they clearly separate at higher values. The uncorrelated distributions descend 
very rapidly while the correlated distributions are order of magnitude larger. 

In Fig. 7 we observe that the full lines and dashed lines, obtained with correlations 
containing only the scalar term, are well separated in the high momentum region. On 
the contrary, the full and dashed-dotted lines, the last ones obtained using correlations 
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without the tensor channels, differ only in a small momentum region around 2.0 fm~^. 
This fact leads us to think that the effects of the tensor correlations are well localized in 
momentum space. 

To investigate better this issue we have performed calculations of the momentum 
distributions using correlations with only the tensor channels active. The result for the 
^^C nucleus are presented in Fig. 8 (dotted line) and is compared with the momentum 
distribution calculated with the full correlation (full line) and with that obtained with 
central and tensor terms only (dashed line). 

We observe that the fiUing of the dip around 2 fm~^ is produced by the sum of the 
effects of the scalar and tensor components. The other central channels are responsible 
of the increase of the tail of the distribution at higher momentum values. 

4 Conclusions 

We have developed a nuclear model to describe density and momentum distributions of 
doubly-closed shell nuclei explicitly considering the short-range correlations. 

This model is based upon the cluster expansion of the CBF theory, but it retains 
only the set of lowest order diagrams which allows for the conservation of the number of 
particles. 

The model has been tested against the results obtained by finite nuclei FHNC calcu- 
lations performed with the same inputs. The agreement between our results and those of 
the more elaborated theory is very good. 
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We have calculated momentum and density distributions of four doubly closed shell 
nuclei with state dependent correlations taken from nuclear matter FHNC calculations. 

We found that the major effects of the correlations in both density and momentum 
distributions are produced by the scalar part of the correlations. The effect of the correla- 
tion functions used in these calculations is a general lowering of the density distributions 
in the interior region and an increasing of the momentum distribution at high momentum 
values. 

The effect of the state dependent terms of the correlation on the density distribution 
is of opposite sign with respect to the effect of the scalar term only. The results obtained 
with the complete correlation are closer to the uncorrelated results than those obtained 
with the scalar term only. The situation is reversed in the momentum distribution case. 

The tensor correlation term produce effects which seem to be rather localized in mo- 
mentum space. In the momentum distributions, the presence of the tensor correlations is 
noticeable in a small momentum region around the point where the uncorrelated distri- 
butions separate from the correlated ones. 

Before concluding we would like to make two comments about the limitations of the 
model with respect to FHNC calculations. 

A first limitation is of practical type. Contrary to what happens in the FHNC case, 
the numerical effort to perform our calculations grows rapidly with the number of single 
particle states. This makes the calculation of the ^ospi^ momentum distribution extremely 
heavy from the numerical point of view. 

The second limitation has a theoretical aspect common to all the models of the same 
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kind. In the FHNC theory the correlation function and the s.p. basis are related through 
the nuclear hamiltonian by the variational principle. In our model they are two different 
input parameters in principle arbitrary. 

We have performed our calculations using reasonable s.p. bases and reasonable cor- 
relations. The s.p. bases have been taken from literature where their parameters have 
been fixed in order to reproduce some nuclear properties within a mean-field model. This 
could mean that some effect of the short range correlations we would like to describe have 
already been averaged out by this procedure. The correlation function used has been ob- 
tained by a minimization procedure done in nuclear matter. We show in Fig. 3 that finite 
nuclei FHNC minimizations seems to prefer correlation functions with an overshooting in 
the region between 1 and 2 fm, but this is not the characteristic shown by the scalar term 
of the nuclear matter correlation function. 

Because of these theoretical hmitations, we beheve that at the present stage the com- 
parison of the results of our model with experimental data is not very meaningful. The 
validity of our work lies in the evaluation of the relative effects produced by the various 
correlation channels. 

Acknowledgments This work has been partially supported by the agreement CYCIT- 
INFN and the Junta de Andalucia. 
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Appendix 

In this appendix we present the expressions of the equations used to calculate the diagrams 
of Fig.l. 

The functions G of Eqs. (|2^) are defined as: 

Gi = Ui + h){9i + 9^) + Kf2 + U){g2 + g,) + Q{h + h){g, + g^) (27) 



G2 = (f,-f,)(g,-g,) + 3{f2-U){g2-g,) + 6{h-mg,-g,) (28) 

G3 = fm + ^hgA + Qhg, (29) 

+ {f2 + U)[\{g2 + gA)~gi-gz]+Kh + h){g, + g,) (30) 

= (/i - /2 + /s - f\) ig2 + gi) + (/2 + f\) (gi - g2 + gs- g^) 

+ 2(/5 + /6)((75 + ^6) (31) 

GQ = {h-h){g,-gA) + {f2-h){^gA-g,) + ^{h-h)g, (32) 

G7 = (/2 - /4)^3 + (/i- 2/2-/3 + 2/4)^4 + 2(/5-/6)^6 (33) 
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G8 = (f,-f,)(g,-g,) + (5U-h){g2-g,) + mg,-g6) 



(34) 



G9 = f3{g2-g^) + f4{gi-2g2-g3 + 2g^)+2U{g,-g,) (35) 

(36) 

To describe the closed shell nuclei we have used a set of s.p. wave functions of the 
form: 

(Pii^Jr,) = Rii^in) J2 < W2s\jm > Yl,{n,)xs{^)xt{^). (37) 

Mi* 

In the above equation Qi indicates the angular coordinates, / and j the orbital and to- 
tal angular momentum respectively, Rnij{r) the radial part of the wavefunction, Yi^ the 
spherical harmonic and < l^l/2s\jm > the Clebsch-Gordan coefficient. 

In this basis the uncorrelated proton OBDM, the first term of Eq. ([TT|) is : 

pS(ri,r2) = ^Y.Ki,ir,)R:,^{r2)i2j + micosO,,) (38) 

nlj 

where Pi is a the Legendre polynomial and 9i2 is the angle between ri and T2- 

The contribution of the other terms has been calculated using the expression (^) of 
the s.p. wave functions and performing a multipole expansion of each correlation function 
/„ and gn of the G functions: 

/.(ri,r2) = 5:/^(ri,r2)Pi(cosM (39) 

L 

where a = 1, 6. 
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The contribution of the diagrams A,B,C and D has beeen calculated expanding in 
multipole each term composing the expressions from (^) to (pS]). Each of these terms 
contains the product between the correlation functions /„ and gn- In the diagram A we 
found terms of the form: 



E 

h,l3 



2/3 + 1 
2/2+1 



I 

n h h 







Pi.icosew) / dr2rlpl'{r2)fr,irur2)giir[,r2) (40) 

J 



In the diagram B we calculate terms of the form: 



and 



[d'r2pl^{r^,r2)pl'{r2,r[)Un2)gp{rv2) = ^ E (2ji + l)(2j2 + 1) 



<A.l(^lX,^2,.2W) E (2/5 + 1) 



h h h 



13,14,15 







h h h 
\^ J 



rcc 

c/r2r2X\,.,,,(r2)<,z,,,(r2)/,:(ri,r2)^7g(r;,r2) Pjcos^^n') (41) 



E / d'r2pr'''ir^,V2)pr'''ir,,v[)Un2)gpin,2] 



Sl,S2 

1 

47r 



Y: V(2^i + l)(2^2 + l)(2ji + l)(2j2 + lX,z„,.(riX,,,,,(r;) 



E (-i)'^-'^^^^M^^ — 



(2/5 + 1) (2/6 + 1) 



I3,l4,l5,h 
roo 



^^2r2X\,,„,,(r2)i?:f„,„,,(r2)/,°(ri,r2)^7f,(ri,r2) 



J2 jl /s 

V 1/2 -1/2 ; 



/ 



h h h 




/ V 



I2 /e h 




^ ^ / / ^ 

n '6 '3 
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h h h 



h h h 



32 ji 1/2 \ h h h 



In the diagrams C and D the functions describing the correlations depend only on two 
radial coordinates, therefore we can perform the multipole expansion of their product. So 



we shall define 



hai3{r23) = f a{r 23) fp{r 23) 



(42) 



In the diagram C we obtain terms of the form: 



"2 
rco 

/ C^?^3?^3Po(^3)^0 (^2,^-3) -Ph (cos ^11') > 



(43) 



while the diagram D we calculate terms of the form: 



^ / rfV2rfV3P^^^^^(ri,r2)p*(r2,r3)p^^^^^(r3,r;)V(r23) = 

S1,S2 

^ (2ji + l)(2i2 + l)it:^,^,,,(ri)i?^3,,„,,(r;)P^,(cosenO 



E 



"2>'2'J2 "3 



^1 ^2 ^4 







i?(ri,r;) 



(44) 



and 



E / (iV2dV3p^^^^^(ri,r2)p^-^*(r2,r3)p^^^^^(r3,r;)/.«;,(r23) 



Sl,S2,«3 
1 



47r 



^ (2ii + l)(2i2 + lX,,,,,(ri)i?^3,,,,,(ri)PacosenO 



"2''2'J2 "3 
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-1) 



where we have defined: 



h+h+h 



( . 



\ 



32 ji h 

[ 1/2 -1/2 ; 



(45) 



/•CO 



(46) 



21 



References 

[1] J.M. Cavedon et al. Phys. Rev. Lett. 49, 978 (1982). 

[2] E.N.M. Quint et al. Phys. Rev. Lett. 57, 186 (1986); E.N.M. Quint et al. Phys. Rev. 
Lett. 58, 1088 (1987). 

[3] R. Jastrow, Phys. Rev. 98, 1479 (1955). 

[4] V.R.P. Pandharipande and R.B. Wiringa, Rev. Mod. Phys. 51 (1979) 821; 

S. Rosati, From nuclei to particles, Proc. Int. School E. Fermi, course LXXIX, ed. 
A. Molinari (North-Holland, Amsterdam, 1982). 

[5] J.W. Clark, Prog. Part. Nucl. Phys. 2, 89 (1979). 

[6] V.R. Pandharipande, Proceedings of Cargese summer School 1989, Cargese 1989, 
edited by J. Tran Tranh Van and J. Negele (Plenum Press, NY, 1990). 

[7] S.C. Pieper, R.B. Wiringa and V.R. Pandharipande, Phys. Rev. C 46, 1741 (1992). 

[8] G. Co', A. Fabrocini, S. Fantoni and I. E. Lagaris, Nucl. Phys. A549, 439 (1992). 

[9] G. Co', A. Fabrocini and S. Fantoni, Nucl. Phys. A568, 73 (1994). 

[10] F. Arias de Saavedra, G. Co', A. Fabrocini and S. Fantoni, preprint, nucl-th/960413. 

[11] F.C. Khanna, Phys. Rev. Lett. 20, 871 (1968); W.J. Gerace and D.A. Sparrow, Phys. 
Lett. B 30, 71 (1969); C. Ciofi degh Atti and N.M. Kabachnik, Phys. Rev. C 1, 809 
(1971). 

22 



[12] M. Gaudin, J. Gillespie and G. Ripka, Nucl. Phys. A176, 237 (1971). 

[13] O. Bohigas and S. Stringari, Phys. Lett. B, 95 9 (1980); M Dal Ri, S. Stringari and 
O. Bohigas, Nucl. Phys. A 376, 81 (1982); F. Dellagiacoma, G. Orlandini and M. 
Traini, Nucl. Phys. A393, 95 (1983); S. Stringari, M. Traini and O. Bohigas, Nucl. 
Phys. A516, 33 (1990). 

[14] O. Benhar, C, Ciofi degh Atti, S, Liuti and G. Salme, Phys. Lett. B 177, 135 (1986). 

[15] M.V. Stoitsov, A.N. Antonov and S.S. Dimitrova, Phys. Rev. C 47, R455 (1993); 
Phys. Rev. C 48 74 (1993); Z. Phys. A345, 359 (1993). 

[16] G. Co', Nuovo Cimento A108, 623 (1995). 

[17] R.B. Wiringa, V. Fiks and A. Fabrocini, Phys. Rev. C 38, 1010 (1988); A. Fabrocini, 
private communication. 

[18] S.Fantoni and S.Rosati, Nucl. Phys. A328, 478 (1979). 

[19] M.M.Renis, Tesi di Laurea, Universita di Lecce, 1995. 

[20] G. Co' and S. Krewald, Phys. Lett. B127 (1984) 145; G. Co', A.M. Lallena and T.W. 
Donnelly, Nucl. Phys. A469 (1987) 684; J.E. Amaro and A.M. Lallena, Nucl. Phys. 
A537 (1992) 585. 

[21] I.R. Afnan and Y.C. Tang, Phys.Rev. 175 (1968) 1337. 



23 



Tables 





A 


C 


<g(i)> 




^tl,p^t2,p 


< t(1) ■ Tf2)Cfl) > 


0/1 20/9 r, 1 1 




< Q(1)t(1) ■ t(2) > 


(5/1 Tif 25/9 n — 1 ) 


0/1 T-iO/o ■nl ZO/'? ^ — i ) 


<r(l)-T(2)Q(l)T(l)-r(2)> 




5tl,p^t2,p(5 - 4:6t3,p) 






D 


<g(i)> 




^tl,p^t2,pSt3,p 


< t(1) ■ t(2)Q(1) > 




0/1 T-iO/'i ■n( .1^ 0/9 ■n ) 

^ L L ^P^ LO fP \ ^lL^P J 


<Q(1)t(1)-t(2)> 




^tl,p'5t3,p(2 — 5t2,p) 


<r(l)-T(2)g(l)r(l)-r(2)> 


5(5ti,p'5t2,p — 2(5ti,p — 25t2,p 


'5n,p5i3,p(55t2,p - 4) 



Table 1. Isospin traces of eq. ([T8|) for each diagram considered. 
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48 Ca 


P 


59.50 


8.55 


4.36 


0.53 




n 


50.00 


7.74 


4.36 


0.53 



Table 2. Coefficients of the Woods-Saxon potential eq. (pBf ). 
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Figure captions 



Fig.l. Set of diagrams considered in our model. Tlie oriented lines represent the 
uncorrelated OBDM (Eq. p!0|) and the dotted lines represent the correlation functions hn 
(Eq. IID. 



Fig. 2. Set of diagrams considered in our model. This set of diagrams corresponds to 
that of Fig.l. In this figure the dashed lines indicate the dynamical correlations (Eq. 

i- 

Fig. 3. Correlation functions used in the ^^C and '^^Ca calculation of Fig. 4. The 
dashed lines show the gaussian correlations and the full lines the correlations obtained 
with a Euler minimization procedure. 



Fig. 4. Comparison between density (upper panels) and momentum distributions 
(lower panels) calculated with our model (dashed lines) and with a full FHNC calculation 
||1CI|| . We have added a factor 1.2 to the densities calculated with the Euler correlations 
and we have been multiplied the Euler momentum distributions by a factor 10. 



Fig. 5. State dependent correlation terms as a function of the interparticle distance 
r?! . See Eq. (^ for the meaning of the various lines. 
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Fig. 6. Proton density distributions calculated with the correlation function of Fig. 5. 
The dotted lines show the uncorrelated results. The full lines have been obtained with 
the full correlations. The results of the calculations performed without the two tensor 
components of the correlations are shown by the dashed dotted lines, while the dashed 
lines show the results obtained with only the scalar term of the correlation (fi). 

Fig. 7. Proton momentum distributions. The meaning of the various lines is the same 
as in Fig. 5. 

Fig. 8. Proton momentum distribution of ^^C. The full line correspond to the full 
calculation (the same as in Fig. 7). The dotted hne has been obtained using only the 
two tensor components of the correlation and the dashed line with the scalar plus tensor 
components. 
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